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ABSTRACT 

The relativistic plasma flows onto neutron stars that are accreting material 
from stellar companions can be used to probe strong-field gravity as well as the 
physical conditions in the supranuclear-density interiors of neutron stars. Plasma 
inhomogeneities orbiting a few kilometres above the stars are observable as X-ray 
brightness fluctuations on the millisecond dynamical timescale of the flows^ Two 
frequencies in the kilohertz range dominate these fluctuations: the twin kilohertz 
quasi-periodic oscillations (kHz QPOs). Competing models for the origins of these 
oscillations (based on orbital motions) all predict that they should be related to 
the stellar spin frequency*^*', but tests have been difficult because the spins were 
not unambiguously known. Here we report the detection of kHz QPOs from a 
pulsar whose spin frequency is known. Our measurements establish a clear link 
between kHz QPOs and stellar spin, but one not predicted by any current model. 
A new approach to understanding kHz QPOs is now required. We suggest that a 
resonance between the spin and general relativistic orbital and epicyclic frequencies 
could provide the observed relation between QPOs and spin. 



Twin kHz QPOs with frequencies between 300 
and 1,300 Hz occur in more than 20 accreting 
neutron stars'^. In seven of these, during ther- 
monuclear X-ray bursts caused by the explosion of 
material accumulated on the neutron-star surface, 
short-hved osciUations are observed with frequen- 
cies Vburst in the range 270-620 Hz, approximately 
constant but different in each source, that are 
thought to be caused by the neutron-star spin^^. 
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The frequency difference Av between the two kHz 
QPOs is close to either Vburst, or t'6urst/2, depend- 
ing on the source. This observed commensurabil- 
ity has given rise to beat-frequency models^' 
Such models identify the higher- frequency ("up- 
per") kHz QPO with the frequency of orbital 
motion, closely around the neutron star, of the 
plasma at the inner edge of the accretion disk 
(within this radius no stable orbits exist and the 
matter plunges in), and the other ("lower") kHz 
QPO with a rotational beat interaction between 
the upper kHz QPO and the neutron-star spin. 
These models predict /S.v to be equal to the spin 
frequency Vspin- In this interpretation, we have 
either v^urst = Jospin or Vburst = '^Vspin-, depending 
on whether one or two hotspots form on the star's 
surface during the bursts. Testing this requires 
searching for harmonic structure in the burst os- 
cillations, but searches^^ are limited in terms of 
statistics as compared to the case of a true pulsar. 

A problem for beat-frequency models is the 
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fact that Av is neither constant^^''^^, nor exactly 
equal^^'^^ to Vspin (but possible solutions in the 
beat-frequency models do exist^°). This has mo- 
tivated the proposal of several alternative mod- 
els to explain the kHz QPOs''~^°, all of which, 
like beat-frequency models, use orbital motion as 
one of the observed frequencies, but which use 
another mechanism - such as general relativistic 
epicyclic motion, which is very fast in these ex- 
treme gravitational fields to generate the other 
frequency. None of these models predict commen- 
surability between spin and kHz QPO frequencies 
such as beat- frequency models do. It has been 
clear for some time that independent measure- 
ments of the neutron-star spin frequency could 
provide the definitive test of the beat-frequency 
idea. The discovery^^ of the first accreting mil- 
lisecond pulsar in 1998 raised the prospect of such 
a test, but until now no accreting millisecond pul- 
sar has shown the anticipated kHz QPOs. 

On 13 October 2002, a bright new outburst of 
the recurrent X-ray transient and accreting mil- 
lisecond pulsar, SAX J1808.4-3658, started^^. Us- 
ing the Rossi X-ray Timing Explorer satellite, we 
obtained ~700 ks of data covering the outburst 
with daily observations from 15 October to 26 
November 2002. We performed a series of fast 
Fourier transforms to search for kHz QPOs, and 
discovered several. One kHz QPO was monitored 
for 12 days (Fig. 1). Its frequency increased from 
570 to 725 Hz between 15 and 16 October, then 
gradually decreased to a minimum of 280 Hz on 21 
October, suddenly increased again to 400 Hz the 
next day, and then resumed its decrease, to 320 
Hz on 26 October. After this the QPO was not 
detected, presumably owing to the limited statis- 
tics resulting from the drop in count rate (by a 
factor of about 4) by this stage of the outburst. 
On 16 October, a second kHz QPO was detected 
simultaneously at a frequency of A;/ = 195±6Hz 
below the first one (Fig. 1, top panel). In the re- 
maining data, this second QPO was not detected 
with a significance greater than 2.5(T. Additional 
broad noise structures of marginal significance are 
sometimes present around 350, 220 and 120 Hz. 

The frequencies, strengths and coherences of 
these two kHz QPOs in SAX J1808.4-3658, as weU 
as the variations in their parameters as a function 
of source luminosity, are very similar to those of 
the twin kHz QPOs observed in those neutron- 



star low mass binaries for which no pulsations are 
seen in their persistent emission (the non-pulsing 
sources^), so we interpret this as the same phe- 
nomenon. The '-^195-Hz frequency difference be- 
tween the two QPOs is the lowest known (in other 
systems Av ranges between 225 and 350 Hz)^. It 
is far below the 401-Hz spin frequency predicted 
in simple beat-frequency models^'^^, yet there is 
a clear commensur ability: Av is consistent with 
I' spin 1'^- This is reminiscent of those non- pulsing 
sources where Au is near I'burst/'^- Hence our re- 
sults support the spin interpretation of the burst 
oscillations in those sources, as well as the suspi- 
cion that Vburst is always near Vspin (and never 
twice that). Moreover, as reported in a com- 
panion papcr^'', in our observations we also de- 
tected four thermonuclear bursts in which we dis- 
covered burst oscillations at the neutron-star spin 
frequency. This discovery further strengthens the 
idea that the burst oscillations are always at the 
stellar spin frequency. 

We now demonstrate that, unless we are 
wrong about the spin frequency of this pul- 
sar, which seems highly unlikely, our observa- 
tions falsify current beat-frequency models, yet 
so pose a severe challenge to all other current 
kHz QPO models. In beat-frequency models^'^'^^, 
the spin-orbit beat interaction occurs at a fre- 
quency Vbeat = n{Vorbit- Vspin), whcrC 

the orbital frequency, and the positive integer 
n is a symmetry factor accounting for multiple 
symmetrically located spin-orbit interaction sites 
(for example, n = 2 for two magnetic poles). 
Simple beat-frequency models have n = 1 and 
the two kHz QPOs at Vbeat and UorUt, so that 

Au = VorUt- J^beat = ^spin- As UOtcd, this is 

clearly inconsistent with our observations. In gen- 
eral, n is the least common multiple (LCM) of 
the number of stellar and the number of orbital 
interaction sites; n = hCM{nspin, Tiorbit)- What- 
ever the value of these numbers, and even allowing 
the observed orbital QPO to be at either Vorbit or 
norbitVorbit, there is no configuration predicting 
the observed Av = Vspinl'^- 

The only solution consistent with current beat- 
frequency models is that Vgpin is half the observed 
401-Hz pulsar frequency (that is, two hotspots) 
and n = 1 (that is, only one pulsar beam or 
magnetic pole interacts with the orbit). How- 
ever, the very stringent amplitude upper limit 
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of 0.014% root mean square (r.m.s.; E.H.M. et 
al., manuscript in preparation) on any coherent 
signal at 200.5 Hz (only 0.38% of the signal at 
401 Hz) makes this solution very unlikely: a fine- 
tuned special geometry would be required to hide 
the true spin frequency that well. We note that 
the sonic point beat-frequency model^ predicts a 
spectrum of many additional weaker QPO peaks 
(higher-order sidebands), one of which (for exam- 
ple, ahs^h'upper ^ 5i>spin) Or, a particularly close 
match, abs(3i';oiuer ~ '^'^spin)^ where Vupper and 
Hower are the upper and lower kHz QPO fre- 
quency, respectively) could be picked out and iden- 
tified with a QPO that we observe. However, this 
would not explain why, contrary to predictions, 
this peak is stronger rather than much weaker 
than the main peaks. In our view, the coinci- 
dence between Az/ and Vhurstl'^ in several other 
sources strongly suggests instead that the relation 
between kHz QPOs and neutron star spin is one 
where (in SAX J1808. 4-3658 and in those several 
other sources) the peaks are separated by i^spj„/2. 
This has not been predicted by any model. 

Inspired by the relativistic resonance models of 
ref. 9 (see also refs 25-27 for related suggestions), 
we note that there are particular radii in the disk 
where the difference between the general relativis- 
tic orbital frequency {y^ and radial frequency (i/^; 
note that — is the periastron precession fre- 
quency) is equal to Ugpin and v spin/ '2.- For a 1.4- 
solar-mass Schwarzschild geometry, these radii are 
~18 and ~23 km, respectively (the best match to 
the observed frequencies in fact occurs at 1.1 solar 
masses, but see below). If a radiatively or mag- 
netically mediated spin-orbit interaction were to 
set up a resonance at one of those radii leading 
to kHz QPOs at Vfj, and Vr, this could explain the 
observed frequency commensurabilities with no in- 
volvement of a beat frequency. Whether Az/ is 
^spin or i^spinl'^ might depend on whether or not 
the relevant radius is within the inner edge of the 
disk; in support of this, we note that all four sys- 
tems where /S.p « Vgpin/2 are those where Pspin ex- 
ceeds 400 Hz, and vice versa for the remaining 
three cases. That the frequencies can shift and 
the observed commensurabilities are often not ex- 
act might be related to a mechanism similar to 
that which makes the resonances occur off the 
naively expected integer ratios in the numerical 
calculations of ref. 28. We note that this combi- 



nation of ideas (adding the spin as an ingredient 
in the relativistic resonance model) can lead to a 
single framework within which the apparently dis- 
crepant frequency commensurabilities observed in 
black holes (which take the form of integer ratios) 
as well as neutron stars can both be understood 
in terms of resonances at privileged radii in an ac- 
cretion disk, namely those radii where the general 
relativistic epicyclic frequencies (and, in the case 
of the neutron-star systems, spin) are commensu- 
rate. 

A surprising third QPO (Fig. 2) near 410 Hz 
was detected on four occasions, just above the 401- 
Hz pulse frequency. Its frequency increased from 
409.3±0.7 to 413.5±0.2 Hz while the frequency 
of the upper kHz QPO increased from 328±4 to 
390±5 Hz. This third kHz QPO is a phenomenon 
not previously seen, and is probably related to the 
pulsating nature of SAX J1808.4 3658. To some 
extent, it resembles the sideband to the lower kHz 
QPO observed^^ in three non-pulsing kHz QPO 
sources: the 410-Hz QPO might be a sideband 
to the pulsation with a similar underlying mech- 
anism. These lower kHz QPO sidebands were 
suggested^^ to be due to Lense-Thirring preces- 
sion, a predicted - but not yet observed - general 
relativistic wobble of the orbital plane (nodal pre- 
cession). The precession frequency is predicted^° 
to be quadratically related to the orbital frequency 
itself, and in the case of both the QPO sidebands 
and the ~410-Hz QPO this is consistent with ob- 
servations. However, as Lense-Thirring precession 
is prograde, a beat between a pulsar beam and 
a precessing orbit would produce a sideband be- 
low, rather than the observed sideband above, the 
pulsar frequency. Another resonance, perhaps at 
the radius where the general relativistic vertical 
epicyclic frequency matches Uspim might be con- 
sidered as an explanation of this phenomenon. 
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Fig. 1. — The power-density spectra obtained for 
SAX J1808.4-3658 during its 2002 outburst. The 
power is given in units of (r.m.s./mean)^/ (Hz x 10""^). 
These spectra show the presence of the pulsations at 
401 Hz, the QPOs, and broadband noise below a few 
hundred hertz. Tiio date the data were obtained is 
given in each panel. All data during those days (for 
the energy range 3-60 keV) were used to create the 
power-density spectra, except for that for 16 October 
for which only the data between 0:53 and 10:00 UTC 
were used and for the energy range 3.3 39 keV. The 
properties of the QPO were obtained by fitting the 
power-density spectra with a constant for the Pois- 
son level (clearly visible at frequencies > 1,000 Hz), 
and three lorentzian functions to fit the continuum 
at low frequencies, the pulsation frequency at 401 
Hz, and the upper kHz QPO. On 16 October, two 
kHz QPOs are visible with frequencies of 499±4 Hz 
and 694±4 Hz, fractional r.m.s. amplitudes (3.3-39 
keV) of 5.2±0.5% (5ct) and 9.1±0.5% (9cr) and widths 
(FWHM) of 35±12 Hz and 82±15 Hz, which we iden- 
tify as the lower and the upper kHz QPO, respectively. 
In the remaining observations one kHz QPO is de- 
tected which we identify as the upper kHz QPO based 
on its factional r.m.s. amplitude of 8-10%, FWHM of 
about 100 Hz, and smooth relation of its frequency to 
that of QPOs and noise features below 100 Hz. No sys- 
tematic differences in its other properties were found 
when this QPO varied in frequency between 725 Hz 
and 320 Hz. 
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Fig. 2. — The power-density spectrum obtained for 
SAX J1808.4-3658 using the combined data from 18- 
19 October and 22-26 October. The power is given in 
units of (r.m.s./mean)^/(Hz xlO~^). This spectrum is 
for the photon energy range 3-60 keV, and shows the 
presence of the pulsations at 401 Hz, the upper kHz 
QPO, the ~410-Hz QPO, and broadband noise below 
a few hundred hertz. The spectrum was fitted with the 
same function as described in Fig. 1, but with an extra 
lorentzian function to fit the extra QPO at ~410 Hz. 
The extra QPO could be detected significantly on four 
occasions: in the combined 18-19 October data, the 
22 October data, the combined 23-24 October data, 
and the combined 25-26 October data. Combining all 
those data, the QPO was significant at a 6.5a level. 
This QPO had an r.m.s amplitude of ~2.5% (for the 
energy range 3-60 keV) and its width (FWHM) was 
between 2 and 8 Hz. 



6 



